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^ . Abstract. An effective quantum pseudospin-1/2 Hamiltonian for Yb2rM207 {TM — Ti 

and Sn) is obtained in terms of the atomic Kramers ground doublet of the LS coupling and 
the crystalline electric field, which is almost described with = ±1/2. It is calculated 
^ • microscopically as the sum of Anderson's superexchange interaction and the magnetic dipole 

' interaction. It is found that it shows a strong exchange anisotropy. 

a' 

Q I 1. Introduction 

, Geometrical frustration and quantum fluctuation provide a promising route to unconventional 

states in magnetism, including various types of spin liquids characterized by the absence or 
emergence of unconventional orders of chirality or multipoles in three spatial dimensions [H [21 [3] . 
The pyrochlore lattice structure offers a prototypical stage where the geometrical frustration 
^ I prevents the magnetic dipole long-range order (LRO) [Hill [5]. One class of magnetic pyrochlore 

CN ■ oxides [6] which have been most intensively studied is a dipolar spin ice, e.g., Dy2Ti207 

and Ho2Ti207 [6l [71 [H [9]. Recently, quantum variants of the spin ice have attracted great 
interest [iQl [III [12] , which include Tb2Ti207 [111 [HI [l5] , Tb2Sn207 [TS], and Pr2Sn207 ^B^, 
Pr2lr207 [13 [E], Yb2Ti207 [20j, and Er2Ti207 [2Q]. 

In the dipolar spin ice, large (111) Ising magnetic moments (~ WfJ-s) of Dy^^ and Ho^^ 
ions, each of which points either inwards ("in") to or outwards ("out") from the center of the 
- tetrahedron, interact with each other mainly through the magnetic dipole interaction [21]. This 

rS I gives a ferromagnetic nearest-neighbor coupling, forming the "2-in, 2-out" ice rule [22] and 

' leaving macroscopically degenerate ground states because of geometrical frustration. The spin 

dynamics is quite slow and the quantum effects are almost negligible. 

Recently, quantum effects have been observed with inelastic neutron-scattering experiments 
on the spin-ice related compounds, Tb2Ti207 [I3l[ll[l5], Tb2Sn207 [l5], and Pr2Sn207 ^7\. 
Pr2TM207 (TM = Zr, Sn, and Ir) shows no magnetic dipole LRO but a (partial) spin freezing 
at Tf ~ 0.1-0.3 K [m [T71 [181 [IS [23]. Pr2lr207 shows a metamagnetic transition at low 
temperatures only when the magnetic field is applied in the [111] direction |19j . as in the spin 
ice. A chiral spin liquid [3] has been detected through the anomalous Hall effect |24j at zero 
magnetic field without magnetic dipole LRO in Pr2lr207 [TU]. Unlike the dipolar spin ice, the 
Curie- Weiss temperature Tew is antiferromagnetic for the zirconate [23j and iridate |18) . The 
stannate shows a significant level of low-energy short-range spin dynamics in the energy range up 
to a few Kelvin [T7], which is absent in the classical spin ice. Vital roles of the planar components 



o 



have also been experimentally observed in Yb2Ti207 and Er2Ti207 [20j. Obviously, quantum 
fluctuations enrich the otherwise classical properties of the spin ice. 

The mechanism in which the quantum fluctuations indicated by these observations appear 
depend on the rare-earth ions, including the difference in Kramers and non-Kramers ions. A 
relatively small localized magnetic moment reduces the magnetic dipole interaction, allowing 
for an important role of the superexchange interaction, as in the case of Pr^+ ions [111 [T2] and 
also expected for Nd^+, Sm^"*", and Yb^"*" ions because of their small moment amplitudes, 3.3/iB, 
0.7/iS) and 4/2^, respectively, for isolated cases. In the case of Tb^'^ ions, the presence of low- 
energy crystalline electric field levels also plays a role [TO]. Remarkably, the large LS coupling 
and the crystalline electric field for localized 4/ electrons makes the superexchange interaction 
highly anisotropic |1H 112] . In fact, the strength of the anisotropy depends on the properties 
of local magnetic doublets of rare-earth ions, producing the variety of experimentally observed 
low-temperature magnetic properties which are mentioned above. Theoretically, it is expected 
that this anisotropic superexchange interaction drives quantum phase transitions among the 
spin ice, quadrupolar states having nontrivial chirality correlations, and the quantum spin ice 
for non-Kramers magentic doublets. 

As for Kramers doublets, on the other hand, understandings of Yb2Ti207 are rather 
controversial. Yb localized 4/ magnetic moments are believed to almost lie perpendicular 
to the (111) directions, unlike Dy, Ho, and Pr cases. It has a ferromagnetic Curie-Weiss 
temperature ~ 0.65 K and a first-order transition has been observed at 0.24 K with the specific 
heat measurement [25j. It has been reported that a neutron-scattering experiment on the single 
crystal shows a LRO in the low-temperature phase [.26], while a combined study using neutron 
diffraction and Mossbauer and ;uSR spectroscopies and the polarized newtron-scattering [28] 
on the polycrystals have indicated that the magnetic dipole correlation remains dynamic. Recent 
neutron-scattering experiments have clarified a [111] rod scattering intensity in the paramagnetic 
phase [29], which has been phenomenologically analyzed by taking into account the anisotropic 
superexchange interaction [30] • 

In this paper, we present a microscopic derivation of an effective quantum pseudospin-1/2 
model, which takes the sum of the magnetic dipole interaction and the superexchange interaction, 
in terms of atomic Kramers doublets of Yb^"*" ions for the frustrated magnet Yb2TM207 (TM 
= Ti and Sn) on the pyrochlore lattice. 

2. Local Kramers doublets of Yb'^+ ions 

The ground state of an isolted Yb^^" ion is characterized by the eightfold degenerate manifold 
^Ft. In the pyrochlore oxide Yb2Ti207, the D^ii crystalline electric field partially lift the 
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degeneracy, yielding the Kramers ground doublet 

7 1 5 

|cr)z) = -aa\Jz = -a) + /3\Jz = -a) + ja\Jz = --^cr), (1) 

with a = ±, real coefficients a « 0.388, p w 0.889, and 7 w 0.242 [3ll [32], and the 
eigenstate \ Jz = M) of the total angular momentum J^, where the local coordinates (xi,yi,Zi) 



(i = 0, • • • , 3) are site-dependent, as explicitly given in Appendix A In particular, the 
quantization axis Zj is always taken to be parallel to the (111) direction that points from the Yb 
site to the center of one of the two tetrahedrons to which the Yb^"*" ion belongs. This ground 
doublet is energetically well separated from the first excitated doublet by 620 K [31j . 

The ground doublet has the following nonvanishing Ising and transverse components of the 
total angular momentum J 

D{a\r\cT')D = ^J||<<.', (2) 



Figure 1. (a) Local level scheme of O 2p holes and Yb 4/ holes, (b) Vitual hopping processes of 
holes contributing to the superexchange interaction, n (n') and i in the state hjhpk'j represent 
the number of / holes at the Yb site r (r') and that of p holes at the O site. 



D{a\r'y\a')n = ^J±<^„ (3) 

where Jy = 7a^ + — 57^ and J± = 2^/7a'y + 4/3^. On the other hand, it does not contribute 
to the atomic quadrupole moment, as is clear from 

D{a\ { J±, r} \a')D =D {a\J+J+\a')D =d {a\J'J"\a')D = 0. (4) 

Let us asign the Pauli matrices ct^. to a psedospin-i that operates on the subspace for the 
the atomic doublet at the site r. Then, the localized magnetic dipole moment is given by 

irir = g.J^^BJr = ^fJ-B g± {a^Xr + al^Vr) + g\\K^r , (5) 

where g± = gjJ± ~ 4.18 and = gjJ\\ ~ 1-77 are reasonably close to experimental values 4.27 
and 1.79, respectively, for Yb2Ti207 [31j, with gj = 8/7 being the Lande factor. This sharply 
contrasts to the non-Kramers magnetic doublet of Pr^"*" ions, where only contributes to the 
magnetic dipole moment and to the magnetic quadrupole moment |1H 112]. 

3. Magnetic dipole interaction 

The magnetic dipole interaction takes the form 



(r,r') 



rhr ■ rhr' {rhr ■ Ar){Ar ■ rhr' 
3 ~ 



(Z\r)3 (Z\r)5 



(6) 



with Ar = r — r' and the summation J2{r,r') o'^^^ ^,11 the pairs of atomic sites, where the localized 
magnetic moment is given by Eq. ([5]). 

4. Superexchange interaction 

To derive an anisotropic superexchange interaction between Yb magnetic moments, it is 
convenient to take the hole picture, where the local ground state is described by a single 4/ 
hole at each Yb site and no hole at O stites. The local Coulomb repulsion among Yb 4/ holes 
and the energy level A of the O 2p holes compared with that of a Yb 4/ hole provide the largest 
energy scales of the problem, as schematically shown in Fig. [1] (a). In fact, the local ground 
state is modified by the nonlocal electron transfer. The amplitude is described by two Slater- 
Koster parameters Vpf-j^ and Vpf^, which represent the transfer integrals between Px/Vy and 
fx{5z^-r^)/ fy{5z2-r^) orbitals and between and f(^5z2-3r^)z orbitals, respectively [33]. Then, we 
perform the fourth-order strong-coupling perturbation expansion of the f-p electron transfer, by 
fully taking into account the virtual electron transfer processes among O 2p orbitals located at 



the center of the tetrahedron and 4/ orbitals located at adjacent Yb sites, which are shown in 
Fig. [1] (b), in the local coordinate frames. The superexchange Hamiltonian is then obtained as 
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// ~ ^2 $Z 5Z 5Z 5Z ^''TT'i^m[^m2^m2f^r,rni,aifr,m2,cr2f^r\m[,a[fr\m2,a2 

{r,r') mi,m2 m\,m'2 '^l.o'a a[,a'^ 

x{^ + h(RtRr') , ArIRv) , , • (7) 

/\ U ^ / mi,m\^\ui,a2 ^ / rrij^, 7712 0-2 

with V±i = Vpfj^ and Vq = Vpf^j as well as the rotation operator of the angular momentum 
from the global coordinates to the local ones at the Yb site r, which is given in Appendix A[ 



Here, Pr,m,a and fr,m,a represent the creation (annihilation) and annihilation (creation) 
operators of 4/ electron (hole) with the orbital m and the spin a at the Yb site r. R^ is 
the rotation matrix of the total angular momentum from the glocal coordinates to the local, as 



given in Appendix A Projecting Eq. ([7]) onto the space of the doublets, we obtain 



Hex = -Jn.n.Y^[9'^KK'+9^{KK'+^lK>) 
{r.r') 

(a^ (ar' ■ flr^r'^ + [^r ' nr,r'^ C^^') , (8) 

with Jn.n. = ^33(7^*)^ (f s) ^^^'^ dimensionless coupling constants as functions of 

gW = 1 - sVQx - -x^ - sVGx^ + — x^, (9) 
2 16 

= 1 + 4\/6x + - 3\/6x3 + -^x^ (10) 
2 lb 

/ = -2(1- 2\/6x + 9x2 _ 3^^3 _^ 9^4^ ^ ^^^-^ 
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5^ = 2^/2fl + ^/6x-^x2 + ^V6x3-^x4V (12) 
V 4 4 8 / 

Here, we have introduced ar = {a^,ar) as well as two unit vectors firy = {cos (pry , — sin (pr,r') 



and n'^ ^, = (sin 0^^^/, cos <?!)r,r')- choice of the local coordinate frames given in Appendix A 

the phase takes (t)r,r' = — 27r/3, 27r/3, and 0, when the pair of Yb sites r and r' corresponds 
to that oi R + ai and R + aj with (i,j) = {(0, 1), (2, 3)}, {(0, 2), (1, 3)}, and {(0, 3), (1, 2)}, 
respectively, where ao = -f (1,1,1), ai = |(-1,1,1), a2 = f (1,-1,1), and as = |(1,1,-1) 
with the lattice constant a and an fee lattice vector R. 

Previously, this form of the nearest-neighbor Hamiltonian with and without the Jnn. term 
is found from the symmetry analysis for the Kramers doublets [12j and microscopically for the 
non-Kramers doublets [111 [T2] , respectively. The present microscopic calculation of the exchange 
coupling constants for the Kramers doublets of Yb^"*" ions remarkably reveals that the symmetric 
exchange part is ferroic in terms of pseudospins, i.e., "in" or "out", and antiferromagnetic in 
terms of real magnetic moments, unlike the case of Pr moments |lll I12| . and that the four 
coupling constants are of the same order in the magnitude with their ratios controlled by a 
single parameter VpfT^/Vpfcj. 
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Figure 2. Dimensionless coupling constants as functions oiVpfT^/Vpfcj. 
5. Summary 

We have obtained a highly anisotropic effective pseudospin-1/2 Hamiltonian for interacting Yb 
Kramers doublets of Yb2rM207, which includes four adjustable parameters; a and 7 describing 
crystal-field doublets, Jn.n. for the magnitude of whole superexchange terms, and Vpf-j^/Vpffj 
specifying the ratios of aniotropic exchange couplings. Theoretical analyses of the model in 
comparison with experiments will be published elsewhere. 
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Appendix A. Local coordinate frames 

The local coordinate frames are chosen as 



a:o = --^(l,l,-2), yo = --L(_i,i,o), ^0 = ^(1,1,1) 



--^(1,-1,2), yi = --L(-l,-l,0), zi = -^(l,-l. 



-1) 



X2 



1 

71 



(-1,1,2), y2 



X3 



(-1,-1,-2), 



-i=(l,l,0), Z2 = -^(-1,1,-1), 
-^(1,-1,0), Z3 = -L(-1,-1,1), 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



for the Yb sites at R + ai (i = 0, • • • , 3). They are invariant under the twofold rotations about 
X, y, and z axes that include the center of the tetrahedron. In particular, all the local z axes 
given above point outwards from the center of the tetrahedron. The coordinate frame for the 
spins is always attached to that for the orbital space in each case. The rotation matrix of the 
total angular momentum j = I + s with the orbital I and the spin s of a single electron takes 
the form 



Rr = exp 



exp 



exp 



(A.5) 



with ipo = 7r/4, 'do = arccos (^l/\/3j , (pi = 3tt/4, -di = — vr + arccos 
t?2 = — •TT + arccos ( , 923 = — 37r/4, and i?3 = arccos il/Vs) . 



^2 = -7r/4, 



Appendix B. Single- hole /-electron states 

The Kramers doublet given in Eq. ([1]) can be expressed in terms of single /-hole states |/^, 
with the orbital P and the spin using 
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